Logarithmic Integrals of Airy Functions 



Bernard J. Laurenzi 
Department of Chemistry 
UAlbany, The State University of New York 
1400 Washington Ave., Albany N.Y. 12222 



October 29, 2012 



Abstract 



Integrals arising in the Thomas-Fermi (TF) theory of atomic structure 
and which contain logarithms of the Airy functions have been expressed 
in terms of the incomplete Bell polynomials. In keeping with the spirit 
of TF theory closed forms for these integrals are sought. 

1 Introduction 

In the theory of atomic structure due to Thomas and Fermi [T] integrals arise [5] 
which contain the Airy function Ai(x) and its derivative Ai'(x) [3J. A typical 
example of these integrals being 



where A(x) is given by 



The value of the integral (to 8 figures) obtained by numerical methods is - 
0.26363171. However, it would be useful and within the spirit of Thomas 
Fermi (TF) theory if an exact or at least an analytic expression for this integral 
and others like it could be obtained. With this in mind, this paper presents 
methods which attempt to move towards that goal. In anticipation of the work 
which follows we set 



I 




z = l- A(x) 



and rewrite the integral above as 




Integration by parts gives 



(1 - z)[l + 2(1 - z) ln(l - z)] s(z) dz. 



(1) 



Where we see that x(z) is the inverse function for A with argument 1 — z. The 
latter function can be obtained using the Lagrange expansion formula |3] i.e. 



'(*) = E I ^ 



fc=i 



(2) 



with the coefficients a^. given by 

afe 



where 



4 (x) = [1 - A(3J)] /x 2 



x=0 



Using (2), the integrals on the right hand side of (1) are elementary and the 
integral I can be written as 



; E 

k=l 



a k M3)-V»(fc/2 + 3)] 
fc! (fc + 2)(fc + 4) ' 



where ij){k) is the Psi (digamma) function [5]. 
1.1 Evaluation of the ak 

The derivatives within a k can in turn be given by the Faa di Bruno formula for 
the (k-l) s * derivative of a composite function [B] as expressed in terms of the 
incomplete Bell polynomials Bk, n (%i,X2, ■ ■ - Xk- n +i) [Z] 



fc-i 

a * = E 

p=l 



d p u~ k / 2 
duP 



where 



u (i) (0) 



d l u(x) 



dx l 



The expression for a k can be quickly reduced to terms involving the Gamma 
function T(k) [5] to get 



fe-i 



p=l 
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The derivatives of the function u(x) can easily be obtained since the power series 
for Ai'(x) is well known [5] i.e. 



Ai'(x)/At(0) = 1 - tllimv £ g| sin( | (n + 4) 7r) (3 1 / 3 *)". 



/••s)y-i'( 

n=0 

Using the latter expression we get 



— 3 5 / 8 r(2/3) 2 • 



(0) = 



3 5 / 6 T(2/3) 

Finally the integral is given by 



oi/3 r(^)sin(f [i+l]7r) 
J (i+2) 



r _ 8 [V-(3)-y.(fc/2+3)] g 

1 ~ 3V3 (fc+2)(fc+4) V - 1 



fc=l 

where Sk is defined as 

' 3 3 / 4 r(2/3) \'' fc-1 



fe! 

p=0 



£ q^S fe - ljP (SW(0),^)(0), . . . , GW-0(O)), 

(4) 



and the reduced quantities 2^(0) are given by 

In obtaining the expression in (3,4) we have used the homogeneous scaling prop- 
erties [TU] of the incomplete Bell polynomials i.e. 

a n B kin (x 1 x 2 , . . .) = B k!n (axi,ax 2 ,...), 
a k Bk, n (xi,X2, ■ ■ ■) = B kt n(axia 2 x 2 , ■ ■ • ), 

to rewrite the Bell polynomials in terms of the reduced quantities u^(0). 



1.2 The Rate of Convergence of the I Integral 

The infinite sum representation for the integral / in (3) is a slowly varying 
function of k. For example, the value of / given by a partial sum containing 
the first ten terms is only 91 % of the value obtained by numerical evaluation 
of the integral. Although a complete analysis (or estimate) of the convergence 
behavior of this sum does not appear to be possible, it is likely that (3) is an 
asymptotic series representation of the integral /. Nevertheless, in an attempt 
to increase the rate of convergence i.e. accelerate the assumed convergence of 
the sum, we note that the magnitudes (Figure 1) of the summands Ik 

t _ 8 [V-(3)-i/>(fc/2+3)] a 
1 k — ^TJW (fe+2)(fe+4) 
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Figure 1: Ik vs. k 



| k -oo, 



25 30 



in (3) are small for values of k greater than 10. 

Furthermore, when Sk is plotted (Figure 2) versus k for values of k > 10, 
it can be seen to vary like powers of 1/k i.e. 

S k = a/k + b/k 2 + c/fc 3 , 

where regression analysis yields a = 0.751653834, b = 2.25325549, c = -6.815672901. 

Figure 2: Sk vs. k 
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Using the curve fitted expression for Sk the expression for / becomes 

10 

j - [V>(3)-^(fc/2+3)] c j_ j 

1 3 1 / 3 (fe+2)(fe+4) ° 



L i-eniaiihdcr ■ 



k=l 



where 



(io If 10 

T 8 r„J„ \ " [^(3)-y.(fc/2+3)] , . \^ [V-(3)-V>(fe/2+3)U s 

-(remainder — 31/3 [« S "1 / , fc (fe+2) (fe+4) f "•" U ] a 2 fc 2 (5+2) (5+4) 

I fc=l J I fc=l 

-4-r/,r V" [^(3)-^(fc/2+3)] \ i 
+C<(T 3 ^ fe3(fc+2)(fc+4) J J' 

The sums Ci appearing in (5) can be calculated in terms of known [12] . closed 
forms i.e. 

oo 

_ [-0(3)- V-(fc/2+3)] _ 89 | ll n /o1 

°1 — fc(fc+2)(fc+4) ~~ 576 ~T 6 111 V^i 

fc=l 

oo 

n n — V* [j(3j-VKfc/g±3)] _ _349_ _ 11/-/^ _ 1 i / 9 n 
" 2 — fc 2 (fc+2)(fc+4) — 1152 16 1 ' 1 . ' 1 6 1U V Z ^ 

fe=l 

V 2 - [V-(3)-^(fc/2+3)] _ 
°3 — / y fc 3 (fc+2)(fc+4) — 
fe=l 



2423 , 5tt 2 IItt 4 7r 2 ln 2 (2) 
" 9216 384 2880 96 



+ ^ + ^ + {^ + ||}C(3) + |L J4 (l/2) ; 

where Li^ is the 4 tft- order polylogarithm function [IT] and £(z) is the Riemann 
zeta function. The value of the sums given above are small with magnitudes 
-0.0389893588, -0.0191766714 , -0.0146522682 respectively. Using these values 
Iremainder = -0.004280449344 yielding a value of -0.2637166702 for I. The latter 
value of I having an error of 0.03 %, this procedure is seen to produce the best 
estimate for an "analytic" representation the integral thus far. 

In order to make further progress in obtaining a closed form expression for 
integrals of the kind sought here, closed form sums of the type 

fc-i 

(i) (^/ 2 )p B k-i, p (u^ (0) , (0) , . . . ,rt k - p H0)), (6) 

p=0 

where (a) n is the Pochhammer polynomial symbol and the u^'(0) terms have 
been given above, must be found. Expressions for closed form sums similar 
(6) have been given by Mihoubi [T3] and Brauchart [TO] albeit with simpler 
arguments of the Bell polynomials. 
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